Adiabatic melting of two-component Mott-insulator states 
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We analyze the outcome of a Mott insulator to superfluid transition for a two- component Bose gas 
with two atoms per site in an optical lattice in the limit of slow ramping down the lattice potential. 
This manipulation of the initial Mott insulating state transforms local correlations between hyperfine 
states of atom pairs into multiparticle correlations extending over the whole system. We show how 
to create macroscopic twin Fock states in this way an that, in general, the obtained superfluid states 
are highly depleted even for initial ground Mott insulator states. 
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I. INTRODUCTION 

The experimental realization of degenerate atomic 
gases py have opened up the possibility of engineering 
strongly correlated many-body quantum states. Such 
states are fundamental for the development of quan- 
tum technologies such as those implementing quantum 
information protocols and quantum computation. A 
paramount example is the realization of a Mott insulator 
(MI) 0, Q in the lowest Bloch band of an optical lattice 
which can serve as a quantum register Q . Bose-Einstein 
condensates of alkali atoms have also been proposed as 
systems to realize macroscopic superposition states [Sj] or 
number correlated states (also called twin Fock states) 
which are useful for quantum metrology 0, 0j B ■ 

In we proposed a method for engineering twin Fock 
states via a two-component MI to superfluid (SF) transi- 
tion. The proposed scheme manipulates hyperfine states 
of atom pairs pinned to single lattice sites and decoupled 
from one another in the MI regime [i,^. The MI limit 
has already proven to be useful starting point for quan- 
tum state engineering since the motional dynamics are 
frozen out and the system decouples to a set of single 
site Hamiltonians to good approximation. In the case of 
a two-component system with doubly occupied sites spin- 
changing collisions will occur and can be accurately con- 
trolled with external magnetic or microwave fields p^ . 
As a result the spin-dynamics of atom pairs behaves like 
an effective two-level system undergoing Rabi-like oscil- 
lations. In contrast to the usual Rabi-model for the cou- 
pling of a pair of single-particle states this represents a 
coherent coupling between a pair of two-particle states in- 
duced by collisional interactions. These local two-particle 
correlations are then transformed into multi-particle cor- 
relations extending over the whole system by melting the 
MI into a SF. This quantum melting can be experimen- 
tally implemented by adiabatically ramping down the 
depth of the lattice potential. 

The build up of long-range correlations when a periodic 
potential is dynamically lowered is a highly non trivial 
process that occurs rather quickly in a one-component 
system 0, HI, |Tl[. The quantum melting of a multi- 
component system is still an open problem that is even 



more intriguing due to the different phases that appear 
in the T — Q phase diagram depending on the ratios of 
the interparticle interactions and of the tunneling rates of 
each component (l^ . An analysis of the two-component 
symmetric case with all interparticle interactions equal 
can be found in |13j] . In 0\ we studied the adiabatic melt- 
ing of an optical lattice in a two-component system with 
different values for the inter-species and intra-species in- 
teractions. We showed that one can create twin Fock 
states with high fidelity via melting of an initial MI state 
with two atoms per site in different internal states. We 
inspected possible sources of experimental imperfections 
such as atom loss during melting or depletion caused 
by final finite interactions. The possibility of creating 
macroscopic superposition states via the adiabatic melt- 
ing of on-site entangled MI states in an infinitely con- 
nected lattice was also briefly considered in [9J. An ex- 
amination of the adiabatic melting in an infinitely con- 
nected lattice, including the effect of the dynamical rela- 
tive phase acquired by the state can be found in . In 
this paper we extend the results in [3| and analyze the 
melting of on-site superposition states in a lattice with 
only nearest neighbor hopping. We show that the final 
SF states obtained from such on-site superposition states 
are highly depleted even if the initial state is a ground MI 
state. For the melted states we calculate the momentum 
distribution and the momentum space particle number 
correlations. These quantities are experimentally acces- 
sible either by time of fiight expansion of the cloud or via 
shot noise measurements [isj. 

This paper is organized as follows. In section |TT] we in- 
troduce the two-mode Bose-Hubbard model (BHM) and 
the initial Mott insulating states that we consider. In sec- 
tion [iTT] we describe the adiabatic approximation and the 
no-crossing rule that we use to infer the final superfiuid 
states. In section IIVI we illustrate how one can create 
macroscopic twin Fock states if the inter-species interac- 
tion is greater than the intra-species interaction during 
melting. In |V] we show using the symmetry properties 
of the initial states and a series of approximations, that 
an on-site entangled MI ground state melts into a highly 
depleted SF state in the limit of intra-species interaction 
greater than the inter-species interaction. 
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II. TWO-COMPONENT BHM 

In this section we introduce the two-component BHM 
that describes the dynamics of ultra-cold atoms in an 
optical lattice. We discuss its experimental feasibility 
and analyze the symmetries of the Hamiltonian. Finally 
we introduce the initial MI states under consideration. 



A. Optical lattice setup 

Our focus is on a state-dependent optical lattice which 
is sufficiently deep that two hyperfine states a and b of 
the atoms are be trapped in the lowest Bloch band of 
the lattice. The resulting system is then accurately de- 
scribed 0] by the two-component BHM with the corre- 
sponding Hamiltonian (?i = 1) 

H = - Y.{Jaa\a, + J,b%) + UY,n'in\ 

+Y E "^("^ - 1) + y E "'('^^ - 1)' (1) 

i i 

where di(bi) is the bosonic destruction operator for an 
a(6)-atom localized in lattice site i, nf = a|ai and 
fi^ = blbi, while denotes summation over nearest- 

neighbors. The parameter Ja(b) is the tunnelling ma- 
trix element for atoms in state a(6); Va^b) ^nd U are 
the on-site intra- and inter-species interaction matrix el- 
ements, respectively. A crucial feature of the optical lat- 
tice setup is that the ratio between the tunnelling and 
interaction matrix elements in Eq. ([T]) is experimentally 
controllable via the lattice laser intensities. Further in- 
dependent control of the interaction matrix elements U, 
V can be achieved with Feshbach resonances jl6|, or by 
shifting the a and b atoms away from each other us- 
ing state-dependent lattices All of these methods 
of control can be exploited dynamically to permit differ- 
ent regimes of the Hamiltonian in Eq. ([1]) to be unitarily 
connected. Together with the long decoherence times of 
ultra-cold atoms this dynamical manipulation can often 
be implemented on near-adiabatic timescales. 

Based on the considerable flexibility of the optical lat- 
tice setup we make a number of simplifications in our 
analysis. First, we restrict the couplings to the symmet- 
ric case where Ja — Jb ~ J and Va = Vb — V. Second, 
we consider states of the system | ip) in which the lat- 
tice is commensurately filled with two atoms per site so 
N\i;) = N\^) has N = 2M, where N = + fi\) 

is the total particle number and M is the total number 
of sites. Finally, we confine our attention to a ID lattice 
with periodic boundary conditions (i.e. M-\- \ = 1) in ad- 
dition to the spatial homogeneity of the matrix elements 
in Eq. ([1]). This is not only a good approximation for 
the center of a large system, where trapping and bound- 
ary effects have a negligible influence, but is also directly 



applicable to recently demonstrated ring-shaped optical 
lattices 

B. Symmetries 

The most basic conserved quantity of the two- 
component BHM is the total number of atoms Na = 
nf, Nb — J2i either of the internal states a or 

b, respectively. With the additional assumptions outlined 
a number of other symmetries arise which we exploit in 
this work. The homogeneity and periodic boundary con- 
ditions assumed results in translational invariance. To 
exploit this we introduce quasi-momentum operators for 
each component as 

1 

1 

^ i—i 

where k is the quasi-momentum in the first Brillouin 
zone and is the position of the site i. The quasi- 
momentum states spanning the first Brillouin zone have 
k ^ Sk X {-Km, --,0, ..,K} with Sk = 2iT/Md and d 
being the distance between lattice sites. For M even 
K.m = M/2 and K = M/2 - 1, while for M odd 
Km = K = {M — l)/2. The total momentum of the 
system is then P — '}2k^^^k + where = c\ck 
and — dl.dk, and is a conserved quantity. The total 
momentum generates cyclic shifts of the lattice sites, e.g. 
i i-^ i + 1 ior the site indices of all operators, via the 
unitary S = exp(— 2 IPtt/M). Its action on the bosonic 
creation operators is 

S dl S"^ — 'il+ii 

SblS^ = bl,, 

simultaneously for all sites i. The homogeneity of the 
system also gives mirror reflection symmetry in which 
site indices for operators interchange as « <-> M — In 
this case the parity operator R is the conserved quantity 
and implements the transformation 



simultaneously for all sites i. Finally, the homogeneity 
and symmetry of the matrix elements result in color in- 
version symmetry a b. Like parity there is a hermitian- 
unitary operator C which is the conserved quantity and 
implements the transformation 

CdlC = bl 

simultaneously for all sites i. These three symmetries are 
schematically represented in Fig. [Ua)-(c). 
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tal particle number N where Na = N — n and Ni, = n 
using rt = 0, 1, 2, . . . , iV. We denote as 5„ the completely 
symmetric subspace of the direct sum 7i„ © Hn-u, de- 
fined by 5 I V') = I V-), I V') = 1 1/-) and C I t/-) = I V) for 
IV') e "Hn © HAT-n where n = 0, . . . , N/2. As we shah 
now outline for the specific initial states considered here 
we fix TV = 2M throughout and the only relevant sectors 
are those where n is restricted to be even or N/2. 



C. Initial states and MI ground states 

Our approach to exploit the two-component BHM to 
engineer quantum states is based on using two experi- 
mentally accessible initial states. Both of these states are 
generated in the MI limit where interactions are domi- 
nant, i.e. U ^ J and V ^ J, and precisely with two 
atoms localized on every site. By controlling the spin- 
changing collisions a doubly occupied single-component 
MI formed in some additional hyperfine state can be 
transformed with near unit efficiency to the state (lo| 



M 



(2) 



FIG. 1: Schematic representations of (a) cyclic shifting, (b) 
mirror reflection and (c) color inversion symmetry for a ID 
lattice with periodic boundary conditions. The colors repre- 
sent atoms in two different hyperfine states a and b. 



with exactly the same number of a and b atoms in each 
lattice site. Note that | '^ab) is exclusively contained in 
the symmetric subspace Sm- By applying a further it/2- 
Raman pulse to the state | ^'afc) an on-site entangled MI 
state 



Symmetries are by definition related to the presence of 
degeneracies [l8| or more precisely to a degeneracy which 
remains for all values of the couplings in the Hamilto- 
nian - a so-called permanent degeneracy. By using these 
symmetries of H it follows that any of its eigenstates 
which are connected by a combination of these symme- 
try operators are necessarily degenerate. Moreover, any 
eigenstates which cannot be mapped on to each other by 
symmetry operators will be non-degenerate if the sym- 
metries of the system are exhaustive fl9| . Other types of 
degeneracies can occur when the couplings attain a par- 
ticular value which indicates that at this point the system 
acquires a higher symmetry as compared to other points. 
This enhancement of symmetry frequently occurs in the 
limits of zero or infinite coupling. An important exam- 
ple in this work is when J = where H decouples into a 
set of single-site Hamiltonians and acquires full permuta- 
tional symmetry. The breaking of the full permutational 
symmetry for finite J splits the degenerate eigenstates 
of the J = Hamiltonian as we will see in proceeding 
sections. 

Particle number conservation for each component 
breaks up the system Hilbert space into sectors where 
N^li/;) ^ Na\ Tp) and Nb\Tp) = Nb\il}). The evolution 
of each component of an initial state in each of these 
sectors can therefore be treated independently. We shall 
use Tin to denote the sectors which possess the same to- 



-bb) — 



(3) 



can be created. This state is an equal superposition of 
all 2*^ possible combinations of the atom- pair states 
I aa) and | hh) over M sites. The states in the superposi- 
tion are contained in the symmetric subspaces Sin with 
n = 0, 1, . . . , \_M/2\ , where [-J denotes the integer part. 
Specifically, in each sector Tiin we form the state 

iV'pr)=ip|ni««). n 1% 



n<j<M 



where P(-) stands for the normalized superposition of all 
permutations of the lattice sites. Since the state | V'p") 
possesses full permutational symmetry it is manifestly 
cyclic- and mirror-symmetric. If we then symmetrize the 
color as [331 



1 



ps 



p.> + l»pV'». (4) 



we obtain states which are each contained in Sim re- 
spectively. The state | "i/aa+bb} can then be seen to be a 
binomial superposition of the states | ^Pp") for all even-n 
symmetric subspaces as 



1 



-bb 



[M/2\ 

E 

n=0 
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The two initial states in Eqs. ([2]) and ([3]) are exact sym- 
metric ground states of the two-component BHM in the 
Umit J = for different regimes of U /V . This is the MI 
regime where the particle-hole spectrmxi is gapped. In 
the case U < V where intra-species interaction is dom- 
inant, and in the subspace 5jv/2j the on-site state | ab) 
is energetically favored making | the nondegenerate 
ground state. This is easily confirmed by comparing it 
to all other configurations and moreover establishes that 
the ground state remains nondegenerate for non-zero J. 

For the opposite case U > V where the inter-species 
interaction is dominant "phase" -separated configurations 
composed of on-site states | aa) = \ o) and | &&) = | •) are 
energetically favored. Thus at J = the full permuta- 
tional and color symmetry results in 2(^^) degenerate 
ground states [sij spanned by all configurations of the 
n states | o) and N/2 — n states | •) distributed amongst 
the M = N/2 sites. Some of these states can be mapped 
by the cyclic shift symmetry S, the mirror R or the color 
C transformation. We denote as G2n C S2n the sym- 
metrized (with respect to C, S and R) ground state man- 
ifold of the BHM in Eq. (P) for J = and [/ > F with 
dimension denoted by g{n, M). 

On this basis we see that each of the states | ^'p" ) com- 
posing I ^aa+bb) IS ouc of many possible ground states in 
the symmetric subspace G2n- However, as will be de- 
scribed in detail in Sec. IV Al once J is non-zero this de- 
generacy is partially split resulting in | '^aa+bb) mapping 
to a superposition of numerous low-lying excited states. 
This is schematically shown in fig. O The reason for 
this can be readily seen. For < J/U ^ 1 nearest- 
neighbor hopping will reduce the energy of configurations 
like I • • • o o o) where a and b atoms form contiguous re- 
gions, compared to configurations like | • o • o • o). Since 
I ^^aa+bb) Contains contributions from all configurations 
it necessarily maps to a superposition of numerous low- 
lying excited states of the two-component BHM in Eq. ([1]) 
for each subspace S2n when J is non-zero. This feature 
significantly complicates the analysis of this state even 
under adiabatic conditions. 



III. ADIABATIC QUANTUM MELTING 

Here we explore how adiabatic quantum melting which 
is an accessible coherent process that can be applied to 
the system and performs a highly non-trivial many-body 
unitary transformation js^l can be exploited as a means 
of quantum state engineering. For both initial states 
I ^'ab) and I "^faa+bb) wc take the quantum melting as end- 
ing in the non- interacting SF limit where U = = 0. In 
this limit the two-component BHM possesses a nonde- 
generate ground state 



E 



Nc,.N,\ _ (co)^°(rfo) 



sf 



, ,._x), (6) 

for any sector Ti.j\j^ in which all the a and b atoms Bose- 
condense into their corresponding zero quasi-momentum 




J = U/J- 



U = V = 



FIG. 2: Schematic energy-level structure for the two- 
component BHM with only one varying parameter U/J and 
fixed U/V. The no-crossing rule assures that the system 
follows the instantaneous eigenstates and performs no level 
crossings if the states have all the same symmetry, there is 
only one varying parameter and the evolution is performed 
sufficiently slow. 



state. 

A feasible way to analyze the properties of the final 
melted states is to use the adiabatic approximation [23| 
in which the system follows the instantaneous eigenstates 
of the two-component BHM in Eq. ([T|) for each value V/J 
and U/J during the melting. The first requirement for 
the analysis to be tractable is the lack of level crossings. 
This condition is greatly aided by the no-crossing rule [l9| 
which states that if the energy levels of a quantum sys- 
tem are plotted as a function of only one parameter, the 
curves for levels of the same symmetry do not cross. The 
no-crossing rule can be followed by using the flexibility 
of the optical lattice setup to ensure that the dynami- 
cal melting only varies the ratio V/J whilst keeping U /V 
constant. Thus the analysis can proceed by following the 
evolution of eigenstates in the symmetric subspace only. 
Also, the BHM is a non-integrable system that shows a 
chaotic spectrum and thus presents energy level repulsion 
[23]. The work in [2^ provides numerical evidence that 
for small systems the symmetries we consider are exhaus- 
tive for the BHM (i.e. no crossings are seen). We extend 
this assumption to arbitrary sized systems. The energy 
levels may approach each other forming so-called avoided 
crossings, but for infinitely slow dynamics there will be 
no transitions between eigenstates. For finite times the 
Landau-Zener model gives the diabatic transition prob- 
abilities [2l[ between two levels, and this generalizes to 
Brundobler-Elser rule [l^l for a multistate problem. 

In principle, adiabatic evolution requires that the dy- 
namical parameter V/J changes in an infinitely long 
time. A crude estimate of the time scale tr needed 
for adiabatic meltin g o f an optical lattice is given by 
U > A^Kiax/MJ^ [nl, where Knax/J stands for the 
initial value in the MI regime. Importantly the adia- 
batic condition does not scale with the system size but 
with the density N/M. Numerical computations in 
indicate that the adiabatic time scale is on the order 
of tr ~ •iVmax/J'^ for a finite realistically sized one- 
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dimensional system. These values are attainable in ex- 
perimental setups and make the adiabatic approximation 
a reasonable assumption. These estimates are valid for 
the case when ?7 = when Eq. ^ decouples as two one- 
component BHM's. Finding analogous time scales for the 
two-component system with U =/= is still an open prob- 
lem but it is reasonable to expect that the corresponding 
time scales are also attainable in experimental setups. 



IV. MELTING THE PRODUCT MI : ^fat) 

We consider first the melting of the product MI state 
I "i'ab) ■ Our analysis begins by examining the outcome 
of quantum melting in a small system by using an exact 
numerical calculation. The results found then generalize 
to an arbitrarily large system and give rise to a final state 
whose usefulness for interferometry is analyzed. 



A. Exact numerics for a small system 



(a) (b) 




FIG. 3: (a) The spectrum of the two-component BHM H 
for M — 6 sites in the symmetric subspace Se, i.e. with 
Na ^ Nt = 6, as a function of V/J with V/U = 0.1 
throughout. To the far left (not plotted) where V^/J S> 1 
the state | ^ab) is the non-degenerate ground state. To the 
right where V/J — the non-interacting SF state | ^t^f) is the 
non-degenerate ground state. As expected these two ground 
states are adiabatically connected (solid line). The first and 
second excited states (dotted and dashed line, respectively) 
are also shown and become degenerate in the V = limit, (b) 
The overlap | ( ipga \ ^ ) | between the instantaneous ground 
state I V'gs) and the SF state | "^^f) as a function of V/J. 



Here we provide some numerical evidence for the 
claims made in the previous section by performing an 
exact numerical diagonalization of the two-component 
BHM Hamiltonian H for a system with M — 6 sites. 
For the state | Vfab) we need only consider the projection 
of H in the subspace Sq. We then compute the ground 
state and low-lying excitations for a sequence of linearly 
decreasing interaction strengths V/ J with a fixed ratio 
V/U = 0.1 throughout. In Fig. ^a,) the eigenenergies 
are shown confirming that the ground state in Sq re- 
mains nondegenerate during the melting, and that within 
the symmetric subspace level repulsion ensures that no 
level crossings occur consistent with the no-crossing rule. 
The final SF ground state is then | ^^f) in which the 
zero quasi-momentum mode of each component is equally 
occupied. Thus the quantum melting transforms local 
initial correlations between pairs of atoms into correla- 
tions between delocalized and macroscopically occupied 
modes. For this reason the final state | 'i'^f) is called a 
twin-Fock state Q. In Fig. [Sib) the overlap |( V-gs | '^tf ) | 
between the instantaneous ground state | V'gs) during the 
ramp and the SF twin-Fock | ^If) is seen to converge to 
unity as a V/J 0. 



B. Large systems and interferometry 

We expand here the results that were investigated in 
. By using the fact that | ^Pab) is a nondegenerate MI 
ground state and the no-crossing rule the outcome of adi- 
abatic melting found for a small system can be readily 
generalized to any size system as 



(a) 





i*tf) 




N/2, N/2. 



(7) 



FIG. 4: (a) The melting of the product MI | *a6) in to the 
twin-Fock SF state | ^'tt). (b) The two internal states a and 
b represent the arms of an interferometer with the rotations 
Rl, R2 and the phase-shift (j> induced by appropriate laser 
pulses. Finally a cj> dependent observable is measured at M 
from which S can be deduced. 



for arbitrary A'". Quantum melting of the MI state | '^ab) 
therefore provides a direct means of creating a iV-particle 
twin-Fock state | ^'tf), as depicted in Fig. dfa). 

Twin-Fock states |6[ have been proposed as useful in- 
put states in interferometric experiments which could po- 
tentially allow measurement sensitivities which scale as 
cx This is the Heisenberg limit for interferometry 

and improves on the standard limit obtained from an un- 
correlated source where the sensitivity scales as cx 1/ ^/N 
[25I [26l . [27I . A typical experimental scheme is a Mach- 
Zenhder interferometer (MZI) that consists of an initial 
beam splitter and two interferometric paths that acquire 
a phase difference (p and are recombined in a final beam 
splitter, as shown in Fig. ID^b). The beam splitter op- 
erations could be implemented by rapid resonant 7r/2- 
Raman pulses. Finally one extracts the phase difference 
measuring a phase dependent quantity. 

To analyze this two-component system within the con- 
text of interferometry we will use the Schwinger boson 
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formalism for constructing SU(2) operators 

i 

i 

and we denote J — [Jx, Jy, Jz)- In many cases Jz is the 
phase-dependent measurement made in the interferome- 
ter. However, for the twin-Fock state (J) = so a differ- 
ent quantity is required. One possibihty is the next order 
angular momentum J^. Note that inefficient detection of 
such quantities can reduce the sensitivities back to the 
standard limit [28j and indirect measurement schemes of 
are needed p9l |. 
If we use J| as the phase-dependent observable the 
resulting sensitivity can be calculated using error propa- 
gation theory. This gives Acfi ~ AJ^/\d{J'^)^/d(j)\ where 

{Jz)4, and AJf = {{Jz)<t, — (-^2)0)^^^ are the average 
and the spread respectively. To calculate {■)^ one can 
either rotate the state or the operators [2^. Since the 
action of the MZI is equivalent to the unitary operation 
exp(— i0Jj,), applying it to the operators gives 



(J2)cos2(0)-H(ij)sin2(0) 
- sin(0)cos((?!))({J2, J^}), 



and 



{Jt)c^ = (Jf)cos4(0) + (44)sin4(0) 

+ sin^((/)) cos^ {(j)){{Jx J Z-, J zJx} 

+{JlJl} + {kJzJJz}) 

— cos (0) sin(0)({J„ if} -f {Jz, JzJJz}) 

- Sm^{(l>) COs{(l>){{Jz, J^} + {X, JxJzJx}), 

where {A, B} = AB + BA is the anticommutator. For 
states which are zero eigenvectors of Jz, such as twin- 
Fock states, this yields 



= 



2 _ sin^ 0((J4) - (JJ)2) + cos^ 4>{jJlJx) 



4C0S2 (/)(J2 



The minimum value of this quantity occurs aX (j) — Q and 
the sensitivity reduces to A(/)(0) — ^(J^)^^/^ after using 
the equality {JxJzJx) = {Jx) which holds for such states. 

The sensitivity can be expressed entirely in terms of 
the one-particle density matrices p° = {a\dj) and p\j — 
(blbj) using 



{Jx) = 



^,3 



(8) 



For the initial (uncorrelated) MI state | ^ab), which pos- 
sesses no off-diagonal correlations pfj — p\j — Sij, this 



yields a sensitivity at the standard limit A(j) = 1/V2N. 
For the final SF | ^fsf) , which possesses long-range cor- 
relations pfj = p\j = 1, we see that asymptotically 

the Heisenberg limit — 1/ ^/N^/2 + N is recovered. 
Thus the scaling A(/) oc iV~" changes from a = 1/2 to 
a = 1 during the melting. An inspection of Eq. ([5]) shows 
that one needs non-zero values of the long-range correla- 
tions pij in order to have better scaling than N given by 
the first term in Eq. ([8]). We showed in p that one can 
obtain sensitivities approaching the HL scaling even in 
the presence of depletion if [/ = 0. Recent results in [s^l 
indicate that our conclusions about the effect of depletion 
are also valid for U > 0. Moreover, we also showed in 
that one can obtain HL sensitivities (a ~ 1) if particles 
are lost during the melting process. 



V. MELTING THE ENTANGLED MI : | *aa+6b> 

In this section we analyze the melting of the entan- 
gled MI state I 'i'aa+bb) ■ In previous work [l3] we demon- 
strated that this state may melt into a maximally entan- 
gled state if the lattice is infinitely connected. In a typical 
experimental set-up with only nearest neighbor hopping, 
the initial permutational symmetry for J = is broken 
for finite J and the final SF state is highly excited. We 
begin by again considering the outcome of quantum melt- 
ing in a small system. Finally after applying a number of 
approximations we compute the momentum correlations 
of the final melted state for a ID lattice. 



A. Small system 

For M — 6 the entangled MI state | '^aa+bb) is an equal 
superposition of all 64 possible combinations of the atom- 
pair states I o) and | •) over M sites. Following Eq. ^ it 
can be written as 

l^^aa+bb) - ^[V2\'f%)+Vu\^%)+V3Q\^l,) 



where the phase-separated states in each subspace S2n 
are 



ps/ 
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(9) 
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FIG. 5: (a) The low-energy spectrum of the two-component 
BHM H for M = 6 sites in the symmetric subspace 54 as 
a function of U/J with V/U — 0.1 throughout. The initial 
state j "i/aa+bb) projected in this subspace is | ^Pps)- The adi- 
abatic evolution of this state involves the ground state and 
two lowest-lying excited states. In the SF limit with U — 
the first-excited level is triply degenerate. Consistent with 
the no-crossing rule we see that energy levels within 54 re- 
pel each other and only avoided crossings are seen, (b) The 
momentum distribution (n^) for one component (it does not 
matter which due to color symmetry) of the final adiabatically 
melted state | ^'f) derived from [ ^aa+bb) with M — 6. 



Here we have used S(-) to denote the normalized equal 
superposition of all equivalent configurations under the 
symmetry operations S, R and C. Each phase separated 
state I ^l"^) is a member of G2n and is a superposition 
of all the g inequivalent symmetric states within it (note 
that g(0,6) = 5(1,6) = 1, while g(2,6) = 5(8,6) = 3). 
The weights of this superposition Eq. ([9]) are given by the 
degeneracies due to symmetrization of these inequivalent 
states, which we denote by s{l) where I = 1, .,5. 

We now proceed to melt each of these phase-separated 
contributions separately. For the subspaces Sq and 
1S2 the corresponding states | ^'pg) and | 'i'ps) are non- 
degenerate MI ground states. Analogous to the melting 
of I ^'ab) in the subspace Sg exact numerical diagonaliza- 
tion of the two-component BHM with U > V confirms 
that there are no level crossings and so these two states 
are adiabatically connected to their corresponding sym- 
metric SF ground states 



$1 



'•)) 

,0)) 



0,12\ 
f ' 

2,10\ 



In this notation i— > we omit the dynamical phase e*" ac- 
quired by the state during adiabatic melting [l^l- In 



contrast subspaces G4 and Ge are three-fold degenerate 
and the states | ^pg) and | ^pg) are both composed of 3 
inequivalent MI configurations. This situation becomes 
more typical with increasing system size. In the J = 
limit these configurations are degenerate, however once 
J is non-zero they split in energy. In Fig. ^a,) the ener- 
gies for the ground state and three lowest lying excita- 
tions in 54 are shown for a linear ramping of U/J with 
U/V = 0.1. This, together with a similar result for 
(not shown), confirm that no level crossings occur. For 
this reason we can map the MI configurations via their 
ordering in energy with their corresponding ground state 
and excitations in the SF limit. The lowest energy con- 
figurations adiabatically connected to the symmetric SF 
ground state are 



$(|....oo)) 



1 



$(|...ooo)) ^ \^'f), 



4,8\ 

sf / ; ' 



respectively, while the two remaining configurations map 
to the lowest and second-lowest excitations in the SF 
limit. Quite generally there are many degenerate SF ex- 
citations in the U — V — limit. For the subspace the 
first two lowest-lying excitations are degenerate, whereas 
the first three are for ^4. To determine the appropriate 
state to which each configuration maps to we apply first- 
order perturbation theory (see Appendix [5^ in the inter- 
and intra-species interaction which splits these states and 
identifies the correct excitations. This then gives for Sq 



S(|»»o.oo)) 
$ (I • o • o • o)) 



cx,l/ : 

6 \ 
cx,2/> 



and for Sa we have 



$ (I • • • o • o)) 
S(|»»o..o)) 



cx.l/ ' 
4 \ 



where | 4'ox.i) and | ^'"x 2) denote the lowest and second- 
lowest SF excitations in 5„ (see Appendix [X]). Note that 
this calculation of the lowest excitations works for ar- 
bitrary M, although it becomes cumbersome for larger 
system size as we would need a very large number of SF 
excitations. 

Combining these results we can write the overall final 
SF state which is adiabatically connected to | ^'aa+bb) '-^ 
I *f) with M = 6 as 



1 r 



(1 *.'f + I ^T)) + (1 ^l^'') + I ^'/')) + VEe^-^ (1 ^If) + I ^tf)) + VT2e 



'"4.1 I ^4 



V6e 



6,6\ 



8 



(a) (i^K) (b) (ntnl) 




FIG. 6: (a) Momentum space particle number correlations 
{n^fig) of the final melted state | ^'f) for one component (it 
does not matter which due to color symmetry) for M — 6. 
(b) The inter-species momentum space particle number cor- 
relations (nlAg) for I *f). 

where we have now written explicitly the dynamical 
phases acquired by each of the states in the superposition. 
An observable of experimental significance is the momen- 
tum distribution for either component. In fact the mo- 
mentum distribution (n^ifc) = (^±fc) ^-^t dependent on 
the precise form of the excited states | Vfcx) since all states 
in the same degenerate manifold have the same expecta- 
tion value ((nl , ) = i in this case). The resulting mo- 

^ 3d 

mentum distribution is shown in Fig.O^b), and gives the 
fraction of the population outside the zero-momentum 
mode (i.e. depletion of the Bose-condensate) for either 
component as D = 1/12 8.3%. Both the intra-species, 
{^'±k'^±q), and the inter-species, {ii±k^±q), momentum 
space particle number correlations do depend on the pre- 
cise form of the excited states | ^ex)- Neither the mo- 
mentum distributions nor the momentum space particle 
number correlations depend on the relative dynamical 
phases. In Fig. [6l[a)(b) these correlations are plotted us- 
ing the excited states determined from first-order pertur- 
bation theory with U/J = 0.01 and V/U = 0.1. 

B. Large systems 

The final SF state obtained via adiabatic melting the 
entangled MI state | '^aa+bb) is given by 

1 .1/2 

where | ^'p") i— > | if we assume that U /V is fixed 

during melting and that the symmetries considered are 
exhaustive and thus no level crossings occur. The nu- 
merical calculation for M = 6 in section IV Al hints that 
it is very complicated to obtain the final state for larger 
M which are inaccessible with exact numerics. Resolving 
the exact splitting of the states in each subspace S2n in 
the MI limit and in the SF limit is both complex and 
impractical for an arbitrarily large M . Instead, we intro- 
duce here some approximations that allow us to calculate 



the momentum distributions and momentum correlations 
of the final states. 

In the MI regime for J = and U > V the color- 
symmetric ground state manifold is (*^) degenerate. For 
finite J, the degeneracy is broken and g{n, M) states (ob- 
tained after symmetrization in S, R) are not coupled by 
the evolution. As shown below for large M we can ap- 
proximate the initial MI state | ^'p" ) as an equal super- 
position of the g{n, M) states. In this way, using the adi- 
abatic approximation we can map these g states into the 
lowest g states in the SF regime. We calculate the lowest 
eigenstates of the BHM in the SF regime U = V = and 
their degeneracies f{l) such that f{l) — g. Even with- 
out calculating the exact coefficients of the final g states 
in the SF eigenbasis, we can compute the momentum 
distribution and the momentum space particle number 
correlations of the final melted state to good approxima- 
tion. 

1. MI limit 

If U /V > 1 and J = the phase-separated state | ^'p") 
belongs to the ground state manifold Q2n C S2n- It is 
constructed as an equal superposition of the (^^^) atom- 
pair states I o) and | •) distributed over M = N/2 sites. 
These atom-pair states belong to the ground state man- 
ifold. However, only g(n,M) < (^^^) of them define the 
symmetrized ground-state manifold G2n and are not cou- 
pled by evolution for finite J. We showed in section FV Al 
that e.g. g{3, 6) = 3 and that for M — 6, the state | 'fpg) 
is a superposition of 3 states that evolve into 3 different fi- 
nal SF states. In general, we can calculate the dimension 
of the subspace Q2n, g{n,M), by noticing that it is iso- 
morphic to the number of color symmetric configurations 
in a bracelet (c.f. FiglT]) with n beads of one color and 
M—n beads of the other color. Such mathematical struc- 
tures have been extensively studied [3l|, [l^l and there 
are analytical formulas available for their properties. For 
even M, the number of complementable bracelets with 
M beads and half of them of each color is given by fsill 

.g(M/2,M) = i[F(f ) + 2f -2 + m 

4 

5(M/2,M) = i[F(f ) + 2f -2 + (I)] M even, 

4 

where F(n) = 2k E<i|„ Hn/d){^ti) + H'2n/d)2^-^ and 
<I>(x) is Euler's totient function. For bracelets with odd 
total number of beads M or even M and n 7^ M/2, the 
number of complementable bracelets is given by [32l | 

gin,M) = i[Fold(G,0,{d : d\M n d\n}) + H{n, M)]. 

Here, Fold(G(-, •), x, {ai, . . . , ak}) gives the last term of 
the list {x,G{x,ai),G{G{x,ai),a2), ■ . ■} generated by 
the cumulative application of the function 

\ J, / 
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FIG. 7: (a) Number of MI ground states {(Imi) that belong to 
each of the g{9, 18) = 765 independent states for finite J for a 
lattice with N/2 = iV„ = iV^ = 18 and U >V. (b) Number of 
independent states with a given degeneracy dA//- The total 
number of independent states is equal to g{9, 18) — 765. (c) 
First energy levels (*) in the SF regime {U = V = 0) and their 
degeneracy / (bar plot) for M = 18. The points (o) show the 
quadratic approximation to the dispersion relation. We plot 
only the first states such that f{l) = g for g = 765. 




FIG. 8: (a) The momentum distribution of the final SF state 
I V&f) obtained via adiabatically melting the MI state | "taa+bb) 
with Na ^ Nb = M = 20 = N/2 and U/V > 1. (b) The 
depletion D = 1 - (no) /Na of the mehed state | 'I'f^^) for 
U/V > 1 as a function of the number of particles A^. The 
solid line shows the best logarithmic fit -0.0748(ln N/2f + 
0.8201 In A''/2 — 1.3633 to the calculated values (shown as o). 



to the elements of the input list {ai, . . . , Ofe} and 

{ Af - 2 
( ) if M even, n odd 
I M 
(lm/2j) Otherwise. 

Note that we recover g{2, 6) — g{3, 6) = 3 obtained for 
the states | ^^^) and | ^l^) given in Sec lVAl 

2. Superposition states 

We now need to write our initial state | ^'p" ) as a su- 
perposition of the g{n, M) states that evolve indepen- 
dently for finite J. Using the no-crossing rule we can 
map the initial superposition | at J = into a final 
superposition | ^Pj"). The coefficients of the superposi- 
tion are given by the the number of bracelet configura- 
tions s{l) that can be obtained by mirror inversion R or 
cyclic permutation S of any of the independent states 
I = 1, ..,g{n,M) in the subspace Q2n- We show s{l) for 
M = 18 in Fig. [3 (a). Note that s{l) = (^). Numer- 
ical results in Fig. [7] (b) show that most (77%) of the g 
inequivalent configurations have the same degeneracy s 
and that this behavior increases with size M (for M = 22 
already 92% of the g bracelets present the same degen- 
eracy) . We can thus assume that the number of states 
isomorphic to each independent configuration is constant 
s — (^) /g. Within this approximation, | ") is thus an 
equal superposition of the first g eigenstates of the SF 
Hamiltonian. 



3. SF Ivmit 

In Appendix B we show how to calculate the eigen- 
states I {e;}) of the Hamiltonian ^ in the SF regime 
where U = V = Q. The lowest eigenstates in ascending 
order in energy for large AI can be calculated efficiently 



assuming a quadratic approximation Eap{l) to the non- 
interacting single-particle energy E{1). Due to the color 
symmetry, the SF eigenstates | e/) are degenerate. The 
degeneracies of the states /(/) can be calculated numeri- 
cally. We show in Fig. [7] (c) the eigenenergies of the first 
SF states and their degeneracies for M = 18. 

To obtain the momentum distribution of | ") we do 
not need to compute the exact distribution of the final 
states I 'J'ex.m) in the {| e;)} basis (see the Appendix IX]) 
because all states in the same degenerate manifold have 
the same expectation value (n°) — (n^). We just need to 
calculate the momentum distribution of the final eigen- 
states {|e;)} and their degeneracies /(/) and sum their 
contributions up to the highest excited state such that 

On the other hand, we saw in Sec. IV AI that 
we need the exact form of the final states | ^cx,m) 
in order to obtain the momentum-number correla- 
tions. Note, however, that '}2L=i I ^^^9 I = 
Hin=i (*cx,m I nfcfi, I *cx,m) whcre we sum over a com- 
plete degenerate manifold. Thus, if we make the further 
approximation that the state | ^Pj") fully populates the 
highest excited manifold we can calculate its momentum 
space particle number correlations without knowing the 
exact states | 4'ex,m)- This is justified for large g because 
the exact form of the highest excited states is concealed 
by the lower manifolds that are fully populated. 

We show the momentum distribution of the melted 
state I in Fig. [8] (a). The final state is highly de- 
pleted. A plot of the depletion D ^ 1 — (no) /N of the 
term | Vt^^^) as a function of N is shown in Fig. [5] (b). 
It indicates the behavior of the depletion of the full state 
I ^ff) because this term has the highest binomial coeffi- 
cient in the sum IjlOp and a very large g. One can see that 
the depletion increases very rapidly with the number of 
particles. 

The momentum-space particle number correlations of 
the melted state | 5'f) are shown in figure [H We observe 
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(a) (ni^fi^) (b) Kof,) 




FIG. 9: (a) The momentum correlations (ntAg) for one 
mode for the state | ^tf) obtained via the adiabatic melting 
of 1 *aa+&(,) with = 40 and U/V > 1. (b) The inter-species 
momentum correlations (n^iiq) of | "^i) with A'^ — 40. Note 
that a twin-Fock state (O would only show a central peak in 
both distributions. 



that the population outside the zero momentum mode is 
higher for the inter-species momentum correlations. Our 
calculations clearly indicate that melting two component 
Mott insulators leads to quantum states with a compli- 
cated structure and long range correlations. A number of 
approximations were necessary to arrive at these results. 
We believe that an exact treatment is unlikely to reveal a 
simple structure which allows an intuitive interpretation 
of the melting process. 



VI. SUMMARY 

We have analyzed the quantum melting of correlated 
two-component MI states into the SF regime using the 
adiabatic approximation. Based on the experimental fact 
that the one component MI-SF transition can be per- 
formed near adiabatically, we have applied it to the two- 
mode case. Using the no-crossing rule we have shown 
that quantum melting of a two-component MI in an op- 
tical lattice provides a viable route for engineering twin 
Fock states useful in interferometric experiments. The 
twin Fock states present only macroscopic occupation of 
the final SF ground state although in general this is not 
the case. We have shown that due to the breaking of the 
translational symmetry when there is nearest neighbor 
hopping the final SF connected to a ground initial state 
can be highly depleted even in the limit of adiabatic be- 
havior. 

Using a set of well justified approximations, we have 
calculated the momentum correlations and momentum 
space particle number correlations for the final states ob- 
tained via melting on-site entangled states created in the 
MI regime. Our results clearly indicate that the melted 
SF states present a complicated structure and long range 
correlations. Further work is needed to ascertain if these 
correlations can be exploited for interferometry or other 
tasks. 
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APPENDIX A: SPLITTING DEGENERATE SF 
EXCITATIONS 

In the non-interacting SF limit where U — V — the 
excited states are made up by expelling pairs of atoms 
of opposite momentum from the ground state Eq. 
Of particular relevance to the exact calculation of the 
melting of | "^aa+bb) in a small system are excitations in 
the subspaces Sm, where \Na — Ni,\ = 0, and Sm-2, where 
\Na — Nb\ — A and M — N/2. Here we shall demonstrate 
that the lowest-lying excitations in Sm and Sm~2 are 
split by the inter- and intra-species interaction to first- 
order in V/J and U / J . 

For the subspace Sm the two lowest-lying excitations 
are degenerate. These form an excitation subspace which 
is spanned by states where a single pair of atoms, of either 
component, is expelled from the condensate to equal and 
opposite momenta with the lowest possible magnitude as 

let) = ^(clrf^+cL4)l*sr^'"-^), 

where ± is used to denote the smallest quasi-momenta 
±Sk. The inter- and intra-species interaction matrix ele- 
ments in this excitation subspace are 

" V{M + + U{M + ^) 2U^l ij 

2f/y^r^ V{M + l-^) + UM 

Upon diagonalizing this 2x2 matrix the resulting gap and 
non-degenerate excitations j^'oii) and [^'^2) (which 
will in general be some superposition of the states | ei) 
and I 62)) are obtained in first-order perturbation theory. 

In a similar way for subspace Sm-2 the three lowest- 
lying excitations are degenerate. This excitation sub- 
space is again spanned by states where a single pair of 
atoms is expelled from the condensate to equal and op- 
posite momenta with the lowest possible magnitude as 

I e.) = \ {cidl + ct_4) (1 v,-+M.-3^ ^ I ^M-SM^i^^ ^ 
lea) = ^(cVjvI;M^*-^)+rfL4|*r''^-^)), 
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To first-order in perturbation theory these states are spht according to the diagonahzation of the 3x3 matrix 
I 



UJ2{1 



hi 



6 

TP 



j_ 

M 



6 
IP 



V{M-l + ^) + U{M^^) 



UJ2(1 



j_ 

M 



2 
IP 



M 



2 
A/2 



V{M + l + + U{M - ij) 



resulting in the non-degenerate excitations | ^ ) , 
|*,^l-)and|vl.^^-). 

APPENDIX B: SF STATES 

In the SF regime {U = 0) the ground state in the 
subspace is given by 

|eo> = ^(|*r^^^^") + l*r-^"'^")) (Bl) 

The excited states | e;) are constructed by expeUing pairs 
of atoms of opposite momentum ±fc. We can denote 
them by | e;) = | ng, ■■,nfc) where Uk — ^(n^ + n^) and 
we omit the negative values because = Uk- An effi- 
cient way of generating the lowest eigenstates in ascend- 
ing order in energy for large M is to assume a quadratic 
approximation Eap(k) ~ to the dispersion relation 



E{k) = 2J(1 — cos(27rfc/M)) for the single particle en- 
ergy in a tight-binding lattice. We can construct the 
eigenstates in ascending order of energy as 

K 

no = iV - nfc 

k=l 
K 

Til — m ~ fc^f^fc 

k=2 

Efc— 1 2 
„ q^n„ 

Uk = U,...,L J k = 2,..,K 

where m is an integer that runs from to some maxi- 
mum value Mm- Due to the color symmetry, the excited 
states |e;) show a degeneracy f{l) (see for example the 
first excited degenerate states Eq. (jAip ) which can be 
calculated numerically. 
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